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A simple method for determination and correction of non-linearity in measurement 
systems is presented . The technique is applicable to a wide range of measuring systems . 
The basic concept , an analysis , and a sample application are given . Non-linearity of the 
Goldstone DSS 13 low noise HEMT 2.3-GHz fS-band ) radiometer system results in noise 
temperature measurement errors . These errors are successfully corrected with this method. 


I. Introduction 

The Deep Space Network (DSN) requires a large number of 
measurement systems for engineering and science data. For 
example, microwave Noise Adding Radiometers (NARs) are 
used for system noise temperature measurements [1]* [3] with 
applications for the determination of radio “star” flux density 
for both antenna efficiency calibrations and scientific value. 

Engineering design and test procedures are used to obtain 
linear system measurement performance. A linear system pro- 
vides equal changes in output reading for equal changes at the 
system input, independent of signal level. If a noise diode is 
turned on and off while the microwave radiometer system is 
switched to the antenna (cold sky) and then to the ambient 
load, a linear system will give equal changes in output read- 
ings. A non-linear system does not. 

Linearity of the DSN NAR, a subassembly of the Precision 
Power Monitor Assembly, is checked with this technique. 1 If 


^‘Precision Power Monitor Assembly,” TM 03 115 A, JPL internal docu- 
ment, pp. 1-3 and 1-4, June 1, 1986. 


the measured noise temperature of a small noise diode, while 
the system is sequentially switched to the antenna and then to 
the ambient load, differs by less than 2 percent, the system is 
considered to be sufficiently linear for use as an NAR. 

The purpose of this article is to evaluate an algorithm [4] 
useful for a wide variety of measurement systems to determine 
and correct for this type of error. 

System non-linearities are only a single example of many 
possible measurement error sources. Other sources of error, 
not discussed in this analysis, include such items as incorrectly 
calibrated noise standards and nonrepeatable front end switches 
for microwave radiometer systems. Each separate error source 
should be systematically analyzed to obtain the required over- 
all measurement accuracy. The following analysis provides a 
method for including non-linearity errors in this list. 

II. Theory 

A. Linear Analysis 

For a linear system, assume an input unknown, 7, to be 
determined in terms of constants and an output reading/? 
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7 = A +BR 


0 ) 


where 


7 = system input parameter (for example, tempera- 
ture, K) 

A,B = constants, to be evaluated, defining a particu- 
lar system 

R = system output reading (arbitrary units: K, mW, 
etc.) 

For a microwave receiving system, 

7 = 7 op = system operating noise temperature for a 
receiving system [5] , K 


A typical DSN microwave radiometer system is presently 
configured with a single microwave ambient termination at 
physical temperature T ^ so that Eq. (5) is not applicable. 
Also, Eq. (6) may not be exactly applicable if R x =£0 when 
T x ~ 0. It is possible to simulate T l by terminating the reading 
device input with an ambient termination (at temperature T P o 
and either setting R , = 0 (zero set) or reading R , (bias). In 
this case, T 1 is estimated using T i = T pi /C, where G is the 
overall gain between the microwave system input and the read- 
ing device input termination. Assuming that T e (defined from 
receiver input to the reading device input termination) is known 
from a laboratory measurement, we have, for this case 


B = 


P 4 


+ T e- T 1 


R 4 -*, 


(7) 


= Ti + T e 

7. = source temperature, K 

T e = receiver effective noise temperature, K 


Evaluating the constants for Eq. (1) 


B = 



A 


t 1 R 4 -t 4 r 1 

*4-*l 


or alternatively 


A = T 1 -BR X 



II 

J-H 

l 

Jg 

(8) 


T , •ir <<1 

(9) 


assuming G is much greater than 1 . 


(2) 

For 7j =R X = 0 


(3) 

7 , 7 

*- X ' 

(10) 


o 

ii 

(11) 

(4) 

Eq. (1) is now used to obtain the system temperature on 
the antenna 7 2 , from the system output reading R 2 , usings 


andi? as calculated above. 


where 


^op = T 2 = A + BR 2 


T x , 7 4 = system input calibration parameters 
R j , R 4 = system output readings with inputs T x and 7 4 

For a microwave radiometer system, the solution for B can 
be accomplished by switching the receiver input between two 
matched terminations at known physical temperatures T P 4 
and T pl , and readings R A and R x . The usual assumption is 
that R = 0 when 7=0. For this case 


B. Non-Linear Analysis 

The linear analysis, Eq. (1), will be in error if the sys- 
tem is non-linear. For small non-linearities, assume a quadratic 
solution (see Fig. 1). Other models are also possible; a qua- 
dratic was chosen for this study since it is the next term in a 
power series, it is consistent with previous studies [3] , and 
good results are obtained with existing microwave radiometer 
systems. 


B = 



,4=0 


T = A+BR+CR 2 (12) 

(5) 

There are various practical techniques for measuring A, B, and 
C. For a non-linear system, constants A and B from Eq. (12) 

(6) are not the same as in Eq. (1). Consider the example of a 
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microwave radiometer system with a manually switched noise 
diode which adds a temperature T n to the microwave system 
input noise temperature when turned on. There are five condi- 
tions to consider, with a noise temperature and reading for 
each. Assume that readings/? are taken for each case. Then for 
the five conditions we find the following input values, 7, and 
output readings,/? (see Fig. 1): 

(1) calibrated termination, ND off: T- T X \R =R X 

(2) antenna, ND off: 7= 7 2 ;/?=/? 2 

(3) antenna, ND on: 7= 7 3 = 7 2 + T n \R =/? 3 

(4) calibrated termination, ND off: 7= 7 4 ;/? =/? 4 

(5) calibrated termination, ND on: 7 = 7 5 = 7 4 + T n \ 
R=R S 

Using these conditions with Eq. (12) 


A +BR l +CR\ = T x (13) 

A+BR 2 +CR 2 2 = T 2 (14) 

A+BR 3 +CR\ = T 2 + T n (15) 

a+br 4 + crI = r 4 (16) 

A +BR S +CR] = T 4 + T n (17) 


Solving five equations with knowns 7 J5 7 4 , R x , /? 2 , /? 3 , 
/? 4 , and/? 5 and 5 unknowns, A, B, C, 7 2 and 7^: 


C = 


7-7 

J 4 J 1 




( 18 ) 


B = -CD 


(19) 


.4 = J, -RR, -CR* (20) 


T 2 = A +BR 2 +CR\ 


( 21 ) 


7 4 and T l are determined for the appropriate configuration, 
such as a microwave radiometer system, using the same con- 
cepts illustrated with Eq. (5). 

The A , B , and C constants adjust the non-linear curve to fit 
the equal input temperature level changes, T n = (7 5 - 7 4 ) = 
(7 3 - 7 2 ) for the observed system output readings. This is the 
central concept of the proposed technique. With a curve de- 
fined by the solution of A , /?, and C, the “primary” unknown, 
7 2 , is solved with Eq. (21). 

These results are general, and can be used either to provide 
an estimate of system non-linearity or to apply a correction. It 
is useful to determine the difference between 7 2 provided by the 
linear and non-linear analyses, Eqs. (1) and (12). An assump- 
tion that the non-linear analysis corrects for non-linearity and 
the linear analysis does not can be used to obtain an estimate 
of the non-linearity error. Defining a linearity factor: 

LF = 7 (linear analysis)/7 (non-linear analysis) 

= 7 (computed from Eq. [1] )/7 (computed from 
Eq. [12]) (23) 

Typically, LF is greater than 1 for an NAR system and less 
than 1 for a total power system. The results section has further 
discussion of the linearity factor. 

It is also useful to consider a linearity correction factor 
(CF) for difference measurements based on Eq. (12). This 
is useful for radio source noise temperature measurements. 
Assume that the source measurement corrected with Eq. (12) 
is given by 


7 - T CF 

sc su 


(24) 


where 

CF = B + 2 CR „ + CT (CF is a ratio obtained from 
off su x 

Eq. [12]) 

7 = R -/?«. = uncorrected source measurement, K 

su on off ’ 

R _ = uncorrected off source measurement 

off 

R - uncorrected on source measurement 

on 


T n = A+BR 3 +CR]-T 2 (22) 

where 

n *l-*l-*\ l +R l 

* S - R <- R 3 +R 2 


The best strategy for correction of source noise tempera- 
ture difference measurements depends on the requirements 
and instrumentation. In some cases the correction should be 
made in near real time, during the sequence of measurements. 
In other cases, particularly with small corrections, the correc- 
tion could be made after the sequence of measurements. Con- 
sider an example with B = 1.02, C= 0.0001, 7^ = 1 K, and 
/? off values of 20 K to 40 K with an elevation angle over the 
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measurement sequence. The corrected source noise tempera- 
ture difference value computed with the averaged R off = 30 K is 
T sc = 1.026 K. Individual corrected values would range from 
1.024 K to 1.028 K. The peak “error” in the single correc- 
tion relative to the full range of uncorrected measurements 
is only about 0.01 dB for this example. This indicates that 
a simple single correction may be satisfactory for many 
applications. 

III. Results 

The above equations have been programmed for an IBM 
compatible personal computer (PC) SuperCalc computer pro- 
gram, assuming a single temperature-calibrated ambient micro- 
wave termination and a known receiver noise temperature, T e . 
Noise temperature measurements were made to verify the pro- 
cedures using the Goldstone DSS 13 antenna at 2.3 GHz 
(S-band). A typical data set taken July, 1987, is shown in 
Fig. 2. Data sets 101, the NAR case, and 102, the total power 
case, are time meshed together in order to minimize other than 
nonlinearity effects. For this data, T x « 0, since there is high 
amplifier gain between the antenna and the square law detec- 
tor. The computer program allows finite values for T x and 
R x (“bias”) if desired. 

The NAR case used an IBM compatible PC configuration. 
An S-band HEMT low noise amplifier was followed by a 
mixer, an IF amplifier, a square law diode detector assembly, 
and an analog voltage to digital interface to the PC. The PC 
controlled the periodic on-off switching of the NAR noise 
diode at the HEMT amplifier input. The noise diode required 
for non-linearity calibration was turned on and off manually 
as required. The total power case used a voltmeter connected 
directly to the analog voltage output of the square law detector. 

Eqs. (1) and (12) are used for the linear and non-linear solu- 
tions tabulated in Table 1 . One-sigma statistical measurement 
errors are in parentheses. The A, B, and C constants are quite 
different for the linear and non-linear analyses. This indicates 
significant non-linearity for these configurations. For the non- 
linear analysis system temperature on the antenna T 2 , results 
for both the NAR and total power cases are in good agreement. 

The linear and non-linear analysis difference for the NAR 
case, data set 101 , is shown plotted in Fig. 3 as the non-linearity 
error. The error goes to zero on the ambient load since this is 
a common reference for both the linear and non-linear analy- 
sis. The linearity factor, LF (Eq. [23]), is shown plotted in 
Fig. 4. A ratio of 1.03, obtained on the antenna, indicates 
about 3 percent non-linearity. This linearity factor is included 
in Table 1 for both the NAR and total power cases, indicat- 
ing an approximate equal magnitude for each. Figures 3 and 4 
include the plus and minus 1 -sigma statistical measurement 


error obtained from the 6 independent data measurements. 
Most of the non-linearity of this configuration is due to the 
square law diode detector. This fact was determined by reduc- 
ing the power level to the detector and repeating the data set 
and analysis. Reducing the power level to the square law diode 
reduces the non-linearity error at the expense of increased zero 
drift in the dc amplifiers following the square law diode, and 
less resolution in the analog to digital converters. 

Figure 5 shows a plot of the differences between the linear 
analysis, Eq. (1), and non-linear analysis, Eq. (12), vs. T op for 
both the NAR case, data 101, and total power case, data 102. 
The small differences between the two different cases after 
correction for non-linearity indicate successful removal of the 
non-linear effects between two separate radiometer operating 
modes. 

The zenith system noise temperature of the DSS 13 S-band 
receiving system with a HEMT low noise amplifier was remea- 
sured by station personnel on August 28, 1987, with three dif- 
ferent time-interleaved data types. The data cases were: (1) IBM 
PC NAR with a square law detector and A/D converter; (2) total 
power with a square law detector and digital voltmeter; and 
(3) total power with a digital power meter. The results are 
tabulated in Table 2 and plotted in Fig. 6. This shows poor 
agreement between the data types using the linear analysis, 
Eq. (1). The non-linear analysis, Eq. (12), corrects for system 
non-linearity and provides good agreement between the three 
data types, increasing confidence in the technique. It is as- 
sumed that the difference between the corrected values for the 
July 2 and August 28, 1987, data is partly due to different 
atmospheric conditions and partly to other measurement bias 
errors. 

IV. Conclusion 

A technique for the determination and removal of system 
non-linearity has been analyzed and demonstrated. This in- 
cludes the concept, equations, and results applicable for DSN 
microwave noise temperature measurements. Comparison cor- 
rections for separate radiometer modes of operation resulted 
in very nearly equal system temperature after correction for 
the system non-linearity. The primary source of radiometer 
system non-linearity was shown to be in the square law detector. 
The best remedy for system non-linearity is proper equipment 
design and system power level settings. Unfortunately, this is 
not always achieved. With a given implemented system and 
operating levels, the above techniques provide a useful tool for 
estimating the system instrumentation non-linearity and cor- 
recting for the non-linearity if necessary. A correction may not 
be advisable if it is of the same magnitude as the correction 
measurement statistical error, or if it is small when compared 
to other error sources. 
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Table 1. DSS 13, July 2, 1987, S-band, zenith, noise temperature measurement analysis 
results (1 -sigma accuracy shown in parentheses) 


Radiometer 

Parameter 

Analysis method 

case 

Linear, Eq. (1) 

Non-linear, Eq. (12) 


A 

0 

0 


B 

1.0034 (0.0004) 

1.0413 (0.0010) 


C 

- 

-0.00012 (0.000002) 

NAR (data 
set 101) 

Noise diode, T n , K 
(on antenna) 

56.4 (0.02) 

57.7 (0.04) 


Noise diode, T n , K 
(on ambient load) 

60.4 (0.09) 

57.7 (0.04) 


System noise 
temperature on 
antenna, T 2 , K 

Linearity factor 
on antenna, LF (ratio) 

31.3 (0.02) 

32.3 (0.03) 

1.034 (0.001) 


A 

0 

0 


B 

281.23 (0.54) 

268.36 (0.59) 


C 

- 

11.82 (0.20) 

Total 

power 

Noise diode, T n > K 
(on antenna) 

59.3 (0.07) 

57.7 (0.05) 

with 

digital 

voltmeter 

Noise diode, T n , K 
(on ambient load) 

54.8 (0.04) 

57.7 (0.05) 

(data set 
102) 

System noise 
temperature on 
antenna, T 2 , K 

Linearity factor 
on antenna, LF (ratio) 

33.7 (0.04) 

32.4 (0.05) 

0.959 (0.001) 


Table 2. DSS 13, August 8, 1987, S-band, zenith, noise temperature measurement analysis 
results (1 -sigma accuracy shown in parentheses) 


Radiometer 

Parameter 

Analysis method 

case 

Linear, Eq. (1) 

Non-linear, Eq. (12) 


NAR (data 

System noise 

30.8 (0.02) 

31.9 (0.10) 

set 101) 

temperature on 
antenna, T 2 , K 

Linearity factor 
on antenna, LF (ratio) 


1.034 (0.003) 

Total 

System noise 

32.9 (0.13) 

31.9 (0.19) 

power 

temperature on 



with 

antenna, r 2 , K 



digital 
voltmeter 
(data set 
102) 

Linearity factor 
on antenna, LF (ratio) 


0.969 (0.003) 

Total 

System noise 

31.9 (0.04) 

31.8 (0.10) 

power 

temperature on 



with 

antenna, T 2 , K 



power 

Linearity factor 


0.997 (0.002) 

meter 

on antenna, LF (ratio) 






Fig. 1. Representative diagram of non-linear system input para- 
meter vs. output reading, assuming that T n = (T 5 - T 4 ) = (T 3 -7 2 ) 
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LINEARITY CALIBRATION DATA SHEET (CTS LIN2, 4-30-87) 

OPERATOR C.GOOPSON FREQUENCY, MHz 2295 DATE 7/2/87 

CONFIGURATION 1.8 V DETECTOR: SN 536 

T f , K 0.471 T m {MASER OR LNA), K 10 INTEGRATION TIME, SEC 40 

T N (IF USING NAR), K 76^5 WEATHER CLEAR 

0.00 0.00 0.00 0.00 0.00 0.00 
1. 0,00 0.00 0.00 0.00 0,00 0.00 
R-, (DETECTOR OR POWER METER IN Z q ) 


0.121 0,120 0.120 0.120 0.120 
2. 31.136 31.108 31.119 31.173 31.193 

R 2 (ANTENNA) 


0.333 0,332 0.332 0,330 0.330 0.329 

3- 87.245 87,295 87,275 87.353 87,391 87.478 

R 3 (ANTENNA AND NOISE DIODE) 

1.095 1.093 1.089 1.089 1.084 1.083 

4. 304.887 304.858 304.576 305.718 305.325 305.551 

R 4 (AMBIENT TERMINATION) 

1.291 1.288 1,284 1.284 1,278 1.277 

5. 365.041 365.110 364.961 365.632 365.374 365.916 

R 5 (AMBIENT TERMINATION AND NOISE DIODE) 

6. 22.23 22,40 22.50 22.64 22.77 22.88 

R 6 (VC) 

7. 1734 1740 1745 1753 1759 1806 

(TIME, UT) 

Fig. 2. Typical radiometer noise temperature data set (time meshed NAR and total power 
cases) taken July 2, 1987, at DSS 13, with 2.3-GHz HEMT low noise amplifier 




SYSTEM NOISE TEMPERATURE (T op ), K 


Fig. 3. Plot of the difference noise temperature for the NAR case, 
data set 101, between the non-linear and linear analyses as a 
function of system noise temperature 



Fig. 4. Plot of the system linearity factor, LF (Eq. [23]), for the NAR 
case, data set 101, as a function of system noise temperature 
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DELTA T. 


□ LINEAR ANALYSIS 
O NON-LINEAR ANALYSIS 



SYSTEM NOISE TEMPERATURE (T op }, K 


Fig. 5. Plot of the differences, ^T op , between the NAR and total 
power cases for both linear and non-linear analyses as a function of 
system noise temperature (see Table 1) 



CASE 2 

TOTAL POWER 
WITH 
DIGITAL 
VOLTMETER 


CASE 3 

TOTAL POWER 
WITH 

POWER METER 


Fig. 6. Plot of the August 28, 1987, zenith, 2.3-GHz noise tempera- 
ture measurements using the NAR and total power cases showing 
linear and non-linear analysis results (see Table 2) 




